Probing confined phonon modes by transport through a nanowire double quantum dot 
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Strong radial confinement in semiconductor nanowires leads to modified electronic and phononic 
energy spectra. We analyze the current response to the interplay between quantum confinement 
effects of the electron and phonon systems in a gate-defined double quantum dot in a semiconductor 
nanowire. We show that current spectroscopy of inelastic transitions between the two quantum dots 
can be used as an experimental probe of the confined phonon environment. The resulting discrete 
peak structure in the measurements is explained by theoretical modeling of the confined phonon 
mode spectrum, where the piezoelectric coupling is of crucial importance. 



The electronic confinement of semiconductor 
nanowires provides a quantization of electron motion 
in the radial direction together with a quasi-one- 
dimensional dispersion in the axial wire direction. The 
resulting electronic properties have been intensively 
studied and suggest that nanowires are promising 
candidates for optical, electronic, and thermoelectric 
applications jlHS]. 

In analogy to the strong electronic confinement, the 
phonons are similarly confined in nanowires, leading to 
a variety of modes with a one-dimensional dispersion [4] . 
This has been shown to yield distinct features in the 
excitonic absorption, such as pronounced phonon repli- 
cas and a substantial broadening of the zero-phonon line 
[5l[6]. However, there have been only few transport stud- 
ies where features of electron-phonon coupling could be 
identified. A recent example is the temperature depen- 
dence of the resistance in Si nanowires which follows a 
power law characteristic for phonon scattering [7 . Going 
beyond such averaged properties, we demonstrate in this 
letter that the individual phonon modes resulting from 
the radial confinement of a nanowire can be probed in 
a transport experiment. This allows for a new kind of 
phonon spectroscopy which provides information on the 
energetic location of the phonon modes as well as their 
interaction strength with the electrons. Such information 
is crucial for any modeling of electronic or thermoelectric 
applications. 

In order to probe the individual phonon modes, we 
apply the resonance condition of tunneling between two 
quantum dots (QDs) inside a nanowire [84 i9j. Trans- 
port properties of gate-defined lateral double QDs in two- 
dimensional electron gases have been widely studied (see 
Ref . ITO^ and references cited therein) . The main feature 
is the strong enhancement of transport if two QD lev- 
els are aligned by tuning the corresponding plunger gate 
voltages. The role of phonon scattering in these lateral 
QD systems has been investigated experimentally [lT|[T2] 
as well as theoretically [13], revealing interesting effects 
such as oscillations due to the phonon interference of the 
spatially separated QDs. While these studies essentially 
probed the bulk phonon structure of the substrate, we 
show that the confined phonon modes in the investigated 
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FIG. 1: (Color online) (a) Scanning electron microscope im- 
age of the measured nanowire lying on the gate electrodes 
together with a three-dimensional model of the entire device 
structure, (b) Energy level diagram of the double quantum 
dot for positive bias, (c) Measurement of the few-electron 
double quantum dot stability diagram at a finite bias of +5mV 
and for weak interdot coupling. The baseline of the lower left 
triangle indicates alignment between the two ground state 
levels in the quantum dots G(1,0) and G(0,1). Red arrows 
indicate a change in detuning A and overall energy £ as a 
function of the two plunger gate voltages Vig2 and Vig4- (d) 
Close-up of the resonance edge between G(0,1) and G(1,0) for 
strong coupling (Mgs = —0.55V). The upper-left inset shows 
the measured current /sd , while its derivative with respect to 
V\g4 is displayed in the main panel to highlight the resonance 
lines running parallel to the ground state transition line. 



nanowire appear as satellite peaks to the resonance for 
tunneling between the QDs. 

The experimental setup is shown in Fig. [ija): A cat- 
alytically grown InAs nanowire with a diameter of 50 nm 
is placed on a substrate with a pattern of metallic gates 
(periodicity of 60 nm), covered by an 18 nm SiN film. 
The sample is identical to the one in Ref. ijA^ where the 
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focus was on single QD properties. Here, the first, third, 
and fifth gate (shown in red color) define an electrostatic 
confinement, leading to two spatially separated QDs in 
the wire. The other two gates (shown in green color) 
are used as plungers to tune the QD filling. A bias volt- 
age Vbias = i5 mV is applied across the structure [see 
Fig.[T](b)]. The measurements are performed at 60 mK in 
a dilution refrigerator. The strong confinement in both 
the radial and growth direction allows for a control of 
the number of electrons in each QD down to the last one 
[m [15] . Figure [TJc) shows that varying the potential of 
the two plunger gates generates the typical stability plot 
for double QD systems [10]. 

In Fig.jljd), we focus on the transition region between 
the states G(0,1) and G(1,0), where the energetically low- 
est level of one of the QDs is occupied, while Coulomb 
charging hinders the addition of a second electron to ei- 
ther QD. The bright line, constituting the common base 
of the double triangle, indicates the relation between 
the plunger gate voltages where the levels G(1,0) and 
G(0,1) of the two QDs align, thus enhancing resonant 
tunneling through the double QD as long as the levels 
are located between the electrochemical potentials (jll/r 
in the left /right part of the nanowire, respectively, see 
Fig. [ijb). Below this line, there is also substantial cur- 
rent flow which can be attributed to inelastic processes 
[TT] . Our key observation is the presence of pronounced 
lines running parallel to the ground state resonance line 
in this region. On these lines, the QD levels maintain a 
fixed detuning A from resonance with an average spac- 
ing of (5 ^ 180 — 200 /ieV. These parallel lines are found 
for both bias directions and cannot be attributed to ex- 
cited states which are several meV above the ground state 
levels [15 . In conventional double QD systems based 
on a two-dimensional electron gas, similar features have 
been found [12] and attributed to an interference in the 
electron-phonon interaction with the bulk phonons. 
As discussed below, our findings cannot be explained by 
this interference. Instead they result from the particular 
phonon structure in nanowires. Thus, the observation 
of these current replicas manifests the interplay between 
electron and phonon confinement in the system. 

To model the phonon scattering environment in the 
nanowire, we calculate the one-dimensional phonon 
modes of the nanowire within an isotropic elastic contin- 
uum model [4] . Our InAs nanowires typically have nearly 
perfect wurtzite crystal structure with the wire oriented 
along the (0001) direction. Neighboring stacking layers of 
the wurtzite structure agree with those in a zinc-blende 
structure in (111) direction Since material parame- 
ters for the InAs wurtzite lattice are largely missing, we 
take them from a rotated bulk zinc-blende structure in- 
stead. Calculations of strain in nanowires have shown 
that such a transformation between the two lattices is 
a good approximation [TJ- The only material param- 
eters which enter the calculation of the phonon modes 
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FIG. 2: (Color online) (a) Phonon spectrum of an InAs 
nanowire with radius = 25 nm (compressional modes only). 
The blue/orange lines depict modes which are of axial/radial 
character for ^ ^ 0. (b) Corresponding density of states of 
the phonon spectrum. 

are the mass density pm = 5680 kg/m^ and the longitu- 
dinal {vi = 4410 m/s) and transverse {vt = 2130 m/s) 
sound velocity, which are taken along (111). Free-surface 
boundary conditions are assumed so that the strain van- 
ishes at the surface of the wire. We neglect the influence 
of the substrate and the Ohmic contacts on the phonon 
mode structure. 

Figure |2|a) shows the compressional (dilatational) 
phonon modes [T8] for our nanowire structure as a func- 
tion of the one-dimensional phonon quasi-momentum 
q resulting from the translation invariance along the 
nanowire axis. These modes decouple into purely ra- 
dial and purely axial modes for q ^ 0. The correspond- 
ing density of states (DOS) is shown in Fig. [2|b). The 
strong peaks in the DOS are roughly equidistant with an 
energy separation of about 180 /ieV, fitting the replicas 
in the current measurement well. It should be noted that 
the singularities in the DOS are a unique feature of the 
one-dimensionality of the nanowire. 

Now we focus on the coupling strength between the 
different phonon modes and the electronic states in the 
QDs. This coupling strength can be quantified by the 
phonon spectral density 

= ^ E i^r - Mirsi^ - ^,.) ■ (1) 

Mf^ describes the electron-phonon coupling element of 
the electronic state in dot i to the phonon mode n with 
one-dimensional quasi-momentum q] here, we assume di- 
agonal electron-phonon coupling. J{uj) takes into ac- 
count the effective coupling strength \Ml'^ — Ml'^l be- 
tween the two QDs as well as the phonon DOS addressed 
above. We consider deformation potential and piezoelec- 
tric coupling to the phonons, for details see Ref. 20, The 
electronic states of the QDs are modeled by cylindri- 
cal Gaussians with radial and axial confinement lengths 
ar and a^, respectively, which is appropriate for a har- 
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FIG. 3: (Color online) Phonon spectral density J{oo) (color 
scale in meV) for different scattering mechanisms and param- 
eters: Deformation potential coupling for varying (a) axial 
confinement length and (b) radial confinement length ar. 
Piezoelectric coupling for varying (c) dielectric constant Sg 
outside the wire and (d) radial confinement length ar- The 
corresponding other parameters are fixed to the values dis- 
cussed in the text. 

monic confinement. From the charging spectrum as well 
as previous measurements we choose the values of 
ar = 15.0 nm and = 20.0 nm. Again we use the 
parameters from the In As zinc-blende structure: con- 
duction band deformation potential Dc = —5.08 eV, 
static dielectric constant el = 15.15, and piezoelectric 
constant 614 = -0.115 C/m^ (Ref. [21]). Note that the 
latter value is uncertain as a variety of different values 
are given in the literature. Taking into account the orien- 
tation of the nanowire, this provides the wurtzite values 
^15 = 631 = —633/2 = — ei4/v^ For the calcula- 

tion of the piezoelectric coupling, the boundary condi- 
tion for the piezoelectric field at the surface of the wire 
is of crucial importance. The wires are located on an 18 
nm thick layer of SiN. Below this layer, metallic contacts 
are present, corresponding to £5 ^ 00. In our radially 
symmetric model [20], we approximate these features by 
assuming an outer dielectric material with = 20. How- 
ever, different values are also studied below. We assume 
zero temperature in analogy to the experimental situa- 
tion. Fig. [3] displays the numerical results for J(a;), see 
Eq. g, for different QD sizes and varying values of e^. 

The deformation potential [Fig.|3|a,b)] couples mainly 
to the radial modes via the diagonal elements of the strain 
tensor [20]. Thus only some of the peaks from Fig. [2] lead 
to pronounced electron-phonon coupling, and the exper- 
imentally observed peak spacing cannot be explained by 
restricting to this scattering mechanism. Reducing the 
QD sizes a^,a2, these peaks receive tails at higher ener- 
gies as the matrix elements for larger i.e., higher cj, 
increase with the reduction in QD size. Due to interfer- 
ence between both QDs, the matrix element IM^'^ — Ml'^l 
vanishes if qd/27T is an integer, where d = 120 nm is the 



separation between the QDs. This causes pronounced os- 
cillations in the spectral density J(cj), which have been 
studied before for bulk phonons [121 US]. The corre- 
sponding energy oscillation period depends on the ve- 
locity V of the phonon mode and takes values between 
S = ?iv27T/d ~ 120 /ieV for the acoustic mode for q ^ 
and ^ ^ 70 fieV for the transverse sound velocity. Even 
smaller velocities are possible close to extrema in the dis- 
persion relation. Thus, these energies are too small to 
explain the origin of the experimentally observed current 
oscillations in Fig. [TJd). 

The piezoelectric coupling also couples dominantly to 
the radial modes (mainly via the coefficient 631) for a 
small outer dielectric constant [Fig. |3|c)]. With in- 
creasing e^, the axial modes (main coupling via 615) pro- 
vide distinct peaks at their dispersion extrema, which 
dominate the coupling if becomes of the order of the 
dielectric constant inside the wire In the latter case, 
we recover a spacing of J ~ 180 /ieV, in good agreement 
with the experimental data. A special situation arises for 
mode 1, which is axial but mainly couples via 633 since its 
displacement is approximately constant. Here, increasing 
Eg causes a relative shift of coupling strength to larger 
energies. For increasing radial confinement length a^, 
the coupling strength typically drops for both the defor- 
mation potential and piezoelectric coupling [Fig. [3jb,d)]. 
However, the dependence of the piezoelectric coupling is 
much weaker than that of the deformation potential. 

Figure [4]^a) shows the experimental current as a func- 
tion of the detuning A between the two QD levels for 
the G(0,1)-G(1,0) resonance. Roughly equidistant shoul- 
ders/peaks (with a separation S ^ 180 — 200 fieV) are 
clearly observed in the tail to the right of the main 
peak for both bias directions and strong interdot cou- 
pling {V\gs = —0.55 V and —0.5 V). Upon increasing the 
barrier height {V\gs = —0.65 V), the current peak be- 
comes more symmetric and no side peaks are observed. 

In our theoretical model, we calculate the current 
through the double QD system following Ref. [131 The 
(diagonal) electron-phonon interaction is treated within 
the independent Boson model [22], and the current is cal- 
culated for a (weak) coupling to the left and right con- 
tacts, F/^/^, using a non-perturbative treatment of the in- 
terdot coupling A finite phonon lifetime 7ph = 40 /ieV 
is taken into account following the work of Refs. 23, In 
the following, we assume 1] = 5 /ieV, F^^ = 10 neV, and 
Tl = 90 neV, which fits the measured peak current and 
provides a part of the asymmetry between forward and 
reverse bias (as the model neglects spin, the full asym- 
metry cannot be accounted for). The results are shown 
in the inset of Fig. |4]^a), solid lines, where the main peak 
at A ^ has a pronounced tail at positive A due to 
tunneling combined with phonon emission. The shoulder 
until A ^ 0.25 meV combines the acoustic phonon mode 
(extending to zero frequency) with its first interference 
oscillation peak, the axial mode 2, and the radial mode 
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FIG. 4: (Color online) Current as a function of the de- 
tuning between the two QDs. (a) Experimental data: In 
the weak coupling case (Mgs = —0.65 V, dotted line) a 
single nearly symmetric resonance occurs. For strong cou- 
pling [Vigs = —0.55 V, solid lines, data as Fig. [l]^d), and 
Vig3 = —0.5 V, dashed line] a right tail appears showing os- 
cillatory behavior. For Vigs — —0.55 V both bias directions 
are shown by different colors. The inset shows numerical re- 
sults for our main parameters (solid lines) and for weak cou- 
pling VL = 200 neV, Fl = F^ = 100 /xeV (dotted line), (b) 
Comparison of peak positions indicated by vertical lines: Ex- 
perimental data (solid lines) and theory (dashed lines) for 
different interdot couplings for the G(0,1)-G(1,0) resonance. 
The dotted line shows data for the G(0,2)-G(l,l) resonance. 



3 (these different features are just about discriminable in 
the theory but appear washed out in the experiment). 
For larger detunings, the peaks can be attributed to the 
onset of the higher axial modes with an energy separa- 
tion of ^ ~ 180 /ieV, in agreement with the experiment. 
These features disappear in the weak coupling case with 
Vt <^T (dotted line). In Fig. gb), we show that the 
peak positions agree between experiment and simulation 
for different couplings [24]. The same peak positions are 
also found at other resonances in this double QD [see 
the dotted line in Fig. |4]^b) and further data from the 
same sample presented in Ref. 15 . Thus, these peak po- 
sitions do not depend on the particular electronic states 
involved. This strongly supports our interpretation that 
the observed peak structure is a signature of the phonon 
confinement. 

In conclusion^ we have found a strong interplay be- 
tween the confined electronic and phononic systems in 
catalytically grown nanowires. Our current spectroscopy 
in the few-electron regime of a nanowire-based double 
quantum dot allows us to probe the confined phonon 
environment of the nanowire. The characteristic peak 
structure can be well explained theoretically if the piezo- 



electric coupling mechanism is taken into account. 
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